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The LHC will have unprecedented sensitivity to flavor-changing neutral current (FCNC) top quark 
decays, whose observation would be a clear sign of physics beyond the standard model. Although 
many details of top flavor violation are model dependent, the standard model gauge symmetries 
relate top FCNCs to other processes, which are strongly constrained by existing data. We study 
these constraints in a model independent way, using a low energy effective theory from which the new 
physics is integrated out. We consider the most important operators which contribute to top FCNCs 
and analyze the current constraints on them. We find that the data rule out top FCNCs at a level 
observable at the LHC due to most of the operators comprising left-handed first or second generation 
quark fields, while there remains a substantial window for top decays mediated by operators with 
right-handed charm or up quarks. If FCNC top decays are observed at the LHC, such an analysis 
may help decipher the underlying physics. 

I. INTRODUCTION 

The Large Hadron Collider (LHC) will have unprecedented sensitivity to flavor changing neutral currents (FCNCs) 
involving the top quark, such as t — > cZ. With a it pair production cross section of about 800 pb and after 100fb _1 
of integrated luminosity, the LHC will explore branching ratios down to the 10~ 5 level Q ^j. Flavor changing neutral 
currents are highly suppressed in the standard model (SM), but are expected to be enhanced in many models of 
new physics (NP). Because top FCNCs are clean signals, they are a good place to explore new physics. There are 
important constraints from B physics on what top decays are allowed, and understanding these constraints may help 
decipher such an FCNC signal. In this paper, we calculate the dominant constraints on top FCNCs from low energy 
physics and relate them to the expected LHC reach using a model-independent effective field theory description. 

Flavor physics involving only the first two generations is already highly constrained, but the third generation could 
still be significantly affected. Of course, the new flavor physics could be so suppressed that it will not be observable 
at all at the LHC. However, since the stabilization of the Higgs mass is expected to involve new physics to cancel the 
top loop, it is natural to expect some new flavor structure which may show up in the top quark couplings to other 
standard model fields. Thus, one may expect flavor physics to be related to the electroweak scale, and then flavor 
changing effects involving the top quark are a natural consequence. 

Although there are many models which produce top FCNCs, the low energy constraints are independent of the 
details of these models. The new physics can be integrated out, leaving a handful of operators relevant at the weak 
scale involving only standard model fields. These operators mediate both FCNC top decays and flavor-changing 
transitions involving lighter quarks. Thus, the two can be related without reference to a particular model of new 
physics, provided there is no additional NP contributing to the B sector. The low energy constraints can be applied 
to any model in which top FCNCs are generated and the constraints on the operators may give information on the 
scale at which the physics that generates them should appear. 

Analyses of FCNC top decays have been carried out both in the context of specific models || and using model 
independent approaches However, in most cases the effective Lagrangian analyzed involved the SM fields after 
electroweak symmetry breaking. As we shall see, the scale A at which the operators responsible for top FCNC are 
generated has to be above the scale v of electroweak symmetry breaking. Thus, integrating out the new physics should 
be done before electroweak symmetry breaking, leading to an operator product expansion in v/A. The requirement of 
SU(2)l invariance provides additional structure on the effective operators ||, which helps constrain the expectations 
for top FCNCs. For example, an operator involving the left-handed (t, b) doublet, the SU(2) gauge field, and the right 
handed charm quark, can lead to b — > S7 at one loop, but also directly to a b — > c transition. If we ignored SU (2) l 
invariance, we would only have the b — > sj constraint, and the resulting bound would be different. An important 
feature of our analysis is that, after electroweak symmetry breaking, the resulting operators can modify even SM 
parameters which contribute at tree level to B physics observables, such as \V c b\. 

The organization of this paper is as follows. In Sec. II we introduce the effective Lagrangian relevant for top 
FCNCs. We also explain why some operators can be neglected and introduce conventions used throughout the paper. 
In Sec. Ill we calculate how these operators affect top quark decays and integrate out the W and Z bosons and the top 
quark to match onto the relevant effective theory at the weak scale. In Sec. IV we relate the experimental constraints 



to the Wilson coefficients calculated in Section III, focusing mostly on observables related to B physics. This leads 
directly to predictions for the top branching ratio. Sec. V contains a summary of the results and our conclusions. We 
include an Appendix with details of the calculations. 



II. EFFECTIVE LAGRANGIAN FOR TOP FCNC 



We consider an effective Lagrangian 



(i) 



where the Oi operators involve third and second generation quarks and the 0[ involve the third and first generations. 
Since we are interested in top quark decays, we define Oi and 0[ in the mass basis for the up- type quarks. 
A complete set of dimension-six operators which give a tcZ or icy vertex are 
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The brackets mean contraction of SU(2) indices, Q 3 and Qi are the left-handed SU(2) doublets for the third and 
second generations, t R and cr are the right-handed SU(2) singlets for the top and charm quarks, H is the SM 
Higgs doublet, H = i^H*, and the index a runs over the SU{2) generators. The first lower L or R index on the 
operators denotes the SU (2) representation of the third generation quark field, while the second lower index refers to 
the representation of the first or second generation field. In this basis all of the derivatives act on the Higgs fields. 
We could also consider operators directly involving gluons, but since the indirect constraints on gluonic currents are 
very weak (see, e.g., 0), we restrict our focus to the electroweak operators in Eq. (Q). The form of the operators in 
Eq. (0) after electroweak symmetry breaking are given in the Appendix. 

Throughout the paper we focus on those new operators that contribute to t — > cZ, cry. In any particular model 
there may be additional contributions to Eq. ([!]) that contribute to AF = 1 and AF — 2 processes in the down 
sector (e.g., four-fermion operators). These operators have suppressed contributions to top FCNCs. When we bound 
the coefficients of the operators in Eq. (|J) from B physics, we neglect these other contributions. In any particular 
model these two sets of operators may have related coefficients. Unless there are cancellations between the different 
operators, the bounds will not get significantly weaker. 

There are other dimension-six operators that can mediate FCNC top decays (for example Ir^D" crB^u). But 
these can always be reduced to a linear combination of the operators included in Eq. (^|) plus additional four-fermion 
operators and operators involving Q lQrM H H fields. For instance, operators involving two quark fields and three 
covariant derivatives can be written in terms of operators involving fewer derivatives using the equations of motion. 
Operators involving two quark fields and two covariant derivatives (e.g., Q 3 D^crD^H) can be written in terms of 
operators involving the commutator of derivatives included in Eq. (|2|) plus operators with one derivative and four- 
fermion operators. Finally, operators involving two quark fields and one covariant derivative can be written in a way 
that the derivative acts on the H field, as in Eq. (||), plus four-fermion operators. 

Of the four-fermion operators which appear after the reduction of the operator basis, some are suppressed by 
small Yukawa couplings and can simply be neglected. However, some are not suppressed, and of those, the biggest 
concern would be semileptonic four-fermion operators, like (tc)(l£). These contribute to the same final state as 
t — > cZ — > c£ + £~. (We emphasize Z — * £ + £~, because the LHC is expected to have the best sensitivity in this 
channel [jl], |^|.) However, the invariant mass of the £ + £~ pair coming from a four-fermion operator will have a smooth 
distribution and not peak around mz , so the Z-mediated contribution can be disentangled experimentally. Operators 
with (tc)(jqq) flavor structure also contribute to t — * c£ + £~ or t — > cy at one loop, but their contributions are suppressed 
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by a/(47r). Finally, operators with the Q l<1rH H H structure either renormalize Yukawa couplings, or contribute to 
FCNCs involving the Higgs (e.g., t — > ch), but we do not consider such processes, as explained later. 

Throughout most of this paper we consider each of the operators one at a time and constrain its coefficient. This 
is reasonable as the operators do not mix under renormalization. One exception is that 0^ L and 0\ L mix with one 
another between the scales A and v, so it would be unnatural to treat them independently. Their mixing is given by 
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where ct2 = a/ sm 2 9w is the SU{2) coupling. (The zero in the anomalous dimension matrix is due to the fact that 
custodial SU(2) preserving operator 0\ L cannot mix into the custodial SU(2) violating 0^ L .) So, we will also carry 
out a combined analysis for these two operators. 

We have written the operators in Eq. (^|) in terms of a single SM Higgs doublet. In principle there may be many 
new Higgs scalars, but only those that acquire a vev will contribute to t — > cZ and 07. Since a triplet Higgs vev is 
tightly constrained by electroweak precision tests, we concentrate on the possibility of multiple Higgs doublets. With 
the introduction of extra Higgs doublets, there are more operators of each particular type (0^ L , 0\ L , etc.), one linear 
combination of which gives rise to t — > cZ and 07. There are also several physical Higgs states that can contribute 
in loops in low energy processes. For each type of operator, a different linear combination of couplings enter in low 
energy measurements. However, without cancellations this will only differ from the one Higgs case by a number of 
order one. This allows our results to be applied to the general case of multiple Higgs doublets. 1 Of course, the Higgs 
sector is also relevant to FCNCs involving the Higgs, such as t — > ch, but we do not consider such processes as they 
are more model dependent. 

Once we go beyond models with minimal flavor violation (MFV) JtJ , the possibility of new CP violating phases in 
the NP should be considered. In MFV models, top FCNC is not observable at the LHC. In models such as next-to- 
minimal flavor violation (NMFV) || top FCNCs could be observable and the Wilson coefficients can be complex. It 
is not always the case that the constraints are weaker when the NP Wilson coefficients are real (in the basis where the 
up type Yukawa matrix is real and diagonal). Rather, interference patterns realized in some of the observables mean 
the constraints are weakest when some of the new phases arc different from or it. We shall point out the places 
where phases associated with the new operators can play an important role and how we treat them. 

In addition to the B physics related constraints we will derive in this paper, one can also use constraints from 
electroweak precision observables. However, these bound flavor-diagonal operators strongly, and the flavor non- 
diagonal operators in Eq. (||) which contribute to top FCNCs are far less constrained. For instance, the Of^ L operator 
corrects the W propagator at one loop and so contributes to the T parameter. The loops involve atoic quark, and 
have one insertion of 0\ L and one insertion of Vt s or V c b- Thus, the contribution is suppressed by |Vt s | ~ \V c b\ ~ 0.04 
relative to an insertion of the flavor diagonal equivalent of 0^ L , Q^HlpW Q3. In contrast, when considering low 
energy FCNC processes, 0^ L will be more strongly constrained then its flavor diagonal version. That is, flavor 
diagonal operators are more tightly constrained by electroweak observables than by low energy FCNCs, while the 
off diagonal operators are more tightly constrained by low energy FCNCs. Moreover, the mixing between these two 
classes of operators is small. It occurs at one loop proportional to y\ |V C {,|, where the factor of yb, the bottom Yukawa 
coupling, is due to a GIM mechanism. Thus, we can think of the flavor diagonal and off diagonal operators as 
independent. And so for the purpose of studying top FCNCs, we are justified in neglecting flavor diagonal operators 
and the relatively weak constraints from electroweak precision tests. 



III. WEAK SCALE MATCHING 



In this section we derive how the NP operators modify flavor changing interactions at the electroweak scale and 
derive the effective Hamiltonian in which the t, W , and Z are integrated out. For numerical calculations we use besides 
the Higgs vev, v = 174.1 GeV, and other standard PDG values f§, \V ts \ = 41.0 x 10" 3 (lO) and m t = 171 GeV Q. 



1 One possible exception is if an extended Higgs sector allows Yukawa couplings larger than in the SM, for example, in a two Higgs doublet 
model at large tan j3. Then a Higgs loop may give additional unsupprcsscd contributions when we match to the Wilson coefficients at 
the electroweak scale. 
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A. Top quark decays 



After electroweak symmetry is broken, the operators in Eq. (||) give rise to t — > cZ and t — > 07 FCNC decays. The 
analytic expressions for the partial widths of these decays are given in Eq. (A2) in the Appendix. Numerically, the 
t — ► cZ branching ratio in terms of the Wilson coefficients is 
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The tcj vertex, which has a magnetic dipole structure as required by gauge invariance, is induced only by the left-right 
operators. The branching ratio for t — > 07 is 



B(t-<ry) = ('i^XVxlO- 4 x 



J RL\ 



(5) 



The analogous expressions for f — > u decays are obtained by replacing Q by C| in Eqs. (0) and (||). 

The LHC will have unparalleled sensitivity to such decays. With lOOfb -1 data, the LHC will be sensitive (at 95% 
CL) to branching ratios of 5.5 x 10~ 5 in the t — > cZ channel and 1.2 x 10~ 5 in the t — ► cj channel |fflj. In the SM, 
£>(t — > cZ, C7) are of order a(V c ba m^/m^) 2 ~ 10~ 13 , so an experimental observation would be a clear sign of new 
physics. Equations (|4|) and (g) will allow one to translate the measurements or upper bounds on these branching 
ratios to the scale of the individual operators. 



B. B decays 

Many of the operators in Eq. (|J) modify SM interactions at tree level (this possibility was discussed in ||). After 
electroweak symmetry breaking, 0\ L gives rise to a bWc vertex with the same Dirac structure as the SM, so the 
measured value of V c b (which we denote V^,° b xp ) will be the sum of the two. This allows us to absorb the new physics 
contribution of C l £ L into the known value of V^ xp — in processes where V c b and C£ L enter the same way, the 
dependence on C^ L cannot be disentangled. For example, the SM unitarity condition, V t * b V t d + V* b V c d + V* b V u d = 0, 
would be violated if one simply shifted the SM values by the NP contributions. However, the CKM fits have unitarity 
built in, so the NP contribution to V c b causes a shift in the values of Vt s and Vtd extracted from the CKM fit, V^ 
and Vid ■ Since we cannot measure all CKM elements independently, we have to replace V ts and Vtd by V^f and V^, 
plus modified NP contributions. (Recall that V ts and V t d are only constrained from loop processes where they enter 
together with new physics contributions.) With these redefinitions we can use V° b p , V^ and V^ in the CKM fit, 
and the NP will only have distinguishable effects in SM loop processes. An analogous procedure applies to the t — > u 
contribution to V^ b p , V^f and V^. Some other operators such as C™ R do not generate a bWc vertex with the same 
Dirac structure as the SM. Thus, their contributions to observables from which V c b is extracted may be disentangled 
as discussed in the following. 

At leading order in the Wolfenstein parameter (Cabibbo angle), A, these relations are: 

Vcb = V: b * p + (v 2 /A 2 )CZ L V tb , 
V ub = V:^ + (v 2 /A 2 )C' L \V tb , 

v* = v; s Rt - (v 2 /a 2 ) (cz L v* + c L \v: s ) , 

v; d = v t T - (v 2 /a 2 ) {ci L v; d + c' L \v: d ) . (6) 

The 0™ R (0'l R ) also modifies the bWc (bWu) vertex, but with different Dirac structure from the SM, so its effects 
can be separated from the SM contribution. Finally, G\ L {0'^ L ) gives tree-level FCNC, since it contains a bZs (bZd) 
interaction. 

At the one-loop level, the operators in Eq. (J2J) contribute to b — > s transitions. The constraints from B physics are 
easiest to analyze by matching these operators onto operators containing only the light SM fields at a scale fj, ~ mw- 
We use the standard basis as defined in |12|. Integrating out the top, W, and Z, the most important operators for 



4 



B — > X s -f and £? — > A" s £ + £ which are affected by NP are 

Ogy = [Sy (1 - 75)6] M - 

Oioa = [«y (1 - 75)6] [^75^] • (7) 

For example, the diagram in Fig. [l] gives a contribution from O^x (denoted by ®) to CV 7 . The coefficients of the 
QCD and electroweak penguin operators, O3 1 ... ) io ! are also modified, but their effect on the processes we consider are 
suppressed. 

Summing the relevant diagrams, the contributions of all operators can be expressed in terms of generalized Inami- 
Lim functions, presented in the Appendix. Setting A = 1 TeV, the numerical results are 2 

C 7l {m w ) = -0.193 + (0.810C£ L + 0.179C£ L + 0.310C$ L - 0.236C^ L + 0.004(7^ - 0.003C%) , 

C 9V (m w ) = 2- [1.56 + (-0.562C£ £ + 44.95C£ Z - 0.885C$ L - 1.127(7^ + 0.0466*^ + 0.004C LH )] , 

C 10A {m w ) = £ [-4.41 + (-7.157Q £ - 598tf£ £ + 3.50C% L - 0.004C^)] . (8) 

The first term in each expression is the SM contribution. Note that the 0\ L contribution is large because it is at 
tree level, while O b LRl O^Ri an d O rr are tiny because they are suppressed by m c /mw and so the constraints on 
these will be weaker. In the case of b — + d transitions the NP contribution has to be rescaled by the 0(1/ A) factor, 
\V t * s V ud /V* d V cs \ w 5.6, and d should be replaced with C\. 

C. AF = 2 transitions 

The operators ££ , C ££ , and O rl also contribute to AF = 2 transitions, i.e., neutral meson mixings. Again, 
the contribution from 0\ L is present at tree level, while the other two contribute starting at one-loop order. The 
relevant functions are again listed in the Appendix. The modifications relative to the SM Inami-Lim function can be 
parameterized as So — * So(l + /iMe 2 '" 4 ') for each neutral meson system. Numerically (setting A = 1 TeV), for B®B a a 
mixing, the effect of the t — > c operators is given by 

h Bs e 21 ^ = 800(C£ £ ) 2 + QmC h LL Ct L - 6.84(C£ L ) 2 + 1.55C£ £ - 2.64Q £ - 0.32(C^ L ) 2 - 1.03C£ £ . (9) 

The contributions of the 0\ operators to B®B® mixing is given by replacing Ci with C[ in Eq. (||) and multiplying its 
right-hand side by A. 

The contribution of the Oi operators to B d B d mixing is obtained by multiplying the right-hand side of Eq. (|^) by 
e 1 ' 9 , where (3 is the CKM phase, (3 = a,rg(-V c dV* b /VtdV t * l ). Whereas the contribution of the 0\ operators to B\B d 
mixing is obtained again by replacing Ci with C- in Eq. (g) and multiplying its right-hand side by — e I/3 /A. 

Finally, the 0[ contribution to K°K° mixing is the same as that to B d B d mixing, up to corrections suppressed by 
powers of A. For the Oj c ontribution to K°K° mixing, one has to replace in Eq. (||) each Wilson coefficient Ci by 
Ci + C*e l/3 (see Eq. ( |A20[ ) in the Appendix), and add to it the additional contribution 

A(h K e 2l,TK ) = 2.26 Re(C^ L C^ L ) e l ? - 5.17 \C't L \ 2 e lf} - 8.35 \C% L \ 2 e lf3 . (10) 

These expressions are valid up to corrections suppressed by A 2 or more. 




FIG. 1: A one-loop contribution from 0^ L (denoted by ®) to Oj 1 . 



2 Throughout this paper we will bound Ci(lTeV/A) 2 and quote numerical results setting A = ITeV. 
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IV. EXPERIMENTAL CONSTRAINTS 



In this section we use low energy measurements to constrain the Wilson coefficients of the operators in Eq. (||). 
Throughout we assume that there are no cancellations between the contributions from different operators. 



A. Direct bounds 

The best direct bounds on the operators in Eq. (|J), as summarized in H], come at present from searches for FCNCs 
at the Tevatron, LEP, and HERA. The strongest direct constraints on t — > cZ and t — > uZ come from an OPAL search 
for e + e~ — > tc in LEP II |l3[ ]. The upper limit on the branching ratio Bit — > cZ, ttZ) < 0.137 bounds the LL and RR 
operators. For neutral currents involving a photon, there is a constraint from ZEUS that looked for e ± p — ► e^tX [fj"I|| . 
This bounds B(t — * try) < 0.0059, and is the strongest constraint on the RL and LR operators with an up quark. The 
other bounds come from a CDF search in Tevatron Run I, which bounds Bit — > cy, uj) < 0.032 (ll| and constrains 
the LR and RL involving a charm. We translate these branching ratios into bounds on the Wilson coefficients and 
list them in the first rows of Tables | and y. The LHC reach with 100 fb _1 data, as estimated in the ATLAS study 0] 
is B(t — > cZ, uZ) < 5.5 x 10~ 5 and Bit — > cry, uj) < 1.2 x 10~ 5 . These will improve the current direct constraints on 
the Wilson coefficients by one and a half orders of magnitude, as summarized in the second rows of the tables. 



B. B — ► X„~f and B -> X s l + l~ 

We first consider the constraints from B — ► X s j. At the scale mj, Oj 7 gives the leading contribution. Using the 
NLO SM formulae from Ref. flTql , we obtain 

BiB^ X sl ) = 10~ 4 x (o.07+ 1 1.807 + 0.08H+ 1.81 AC 7l im w )\ 2 ) , (11) 



where ACV 7 (mw) is the NP contribution to Cy 7 at the /i = m^y matching scale. The current experimental average |{L7| , 
B(B A s7 ) = (3.55 ± 0.26) x 10~ 4 , implies at 95% CL 3 (setting A = 1 TeV) 

- 0.07 < Cl L < 0.04 or 1.2 < C£ L < 1.3 , 

-0.3 < C'l L < 0.16 or 5.3 < C£ L < 5.8 , 

-0.2 < C^ L < 0.1 or 3.1 < C% L < 3.4 , 

-0.1 < C b RL < 0.24 or -4.5 < C b RL < -4.1 , (12) 

The first (left) intervals are consistent with the SM, while the second (right) ones require new physics at the 0(1) 
level. The non-SM region away from zero is disfavored by b — > s£ + £~ discussed below, but we include it here for 
completeness. For the operators whose contributions are suppressed by m c , we find 

- 14 < Cl R < 7 , -10 < C b LR < 19 , (13) 

and no meaningful bound for C RR . To obtain the results in Eq. ( |l2|) and (fL3|), we assumed that the NP contributions 
are real relative to the SM, i.e., that there are no new CP violating phases. Had we not made this assumption, the 
allowed regions would be annuli in the complex Cj planes. 

Next we consider B — > X s £ + £~ . The theoretically cleanest bound at present comes from the inclusive B — > X s £ + £~ 
rate measured for 1 GeV 2 < q 2 < 6 GeV 2 @ 

B(B -> X s £+r) 1GcV 2 <q2<6GcV 2 = (1.61 ± 0.51) x lO" 6 . (14) 

Due to the unusual power counting in B — > X s £ + £~ . the full set of C(a s ) corrections are only included in what 
is called NNLL order, achieving an accuracy around 10%. For the SM prediction we use the NNLL calculation as 
implemented in Ref. |19]. This calculation does not normalize the rate to the B — > X£v rate; doing so would not 



3 Hereafter all constraints are quoted at 95% CL, unless otherwise specified. 
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FIG. 2: Constraints from B -» X s -) and B -> X s £+^" in the C£ L - C£ L plane. The red, green, and blue regions denote 68%, 
95%, and 99% CL, respectively. The region between the dashed lines is beyond the LHC sensitivity. 



improve the prediction significantly and would unnecessarily couple different operators' contributions. We include the 
modifications of CV 7 , Cgy, and C\qa due to the new operators at lowest order. With our input parameters, we obtain 

B{B -> X s £ + e-) 1<g2<6GeY2 = 1(T 6 x [l.55 + 35100 [\AC 9V (m w )\ 2 + \AC WA (m w )\ 2 ] + 0.45 \AC 7l (m w )\ 2 

+ Re[(180 + 5i)AC 9y (mvi/)] - 360Re[ACW(mw)] (15) 
-Rc[(0.17 + 0.04i)AC 77 (mvi/)] - 200Re[AC 9V (m w )* AC l7 (m w )] } . 

The simplest way to proceed would be to bound C/ 7 , Cgy , and C\qa separately at fi — mw , assuming that the others 
have their SM values, and use this to constrain new physics. This procedure would not be consistent, since the NP 
necessarily affects these Wilson coefficients in a correlated way. Instead, we directly constrain the coefficients of O l £ L , 
Oll> ®lr-i an d O b LR , which also yields stronger constraints. With A = ITeV, we obtain 

10~ 2 or - 5 x 10~ 3 < C\ L < 3 x 10~ 3 , 
or 1.7<Cg t <3, 

(16) 

The combined constraints from B — » X s ~f and B — > X s £ + £~ on these four Wilson coefficients are shown in Table || 
in the Conclusions. We plot in Fig. || the bound on the LL operators in the C^ L - C\ L plane. The SM corresponds 
to the point (0, 0). A measurement or a bound on the t — > cZ branching ratio corresponds to a nearly vertical band. 
The LHC is sensitive to this whole plane, except for the band between the dashed lines. 

The above bounds were derived assuming that the NP contribution is real relative to the SM. ft is conceivable that 
improved measurements of B — > X s £ + £~ will lead to constraints on the CP violating phases before the LHC is be 
able to probe top FCNCs. Thus we postpone a full analysis with complex NP Wilson coefficients until more data 
become available. 



-1.1 < C U LL <0.3, 
-1.8 x 10~ 2 < C'l L < -1 x 
-0.5 < C% L < 0.7 
-2.0 < C h m < 3.5. 



C. Exclusive and inclusive b — > civ decays 

In this subsection we investigate the constraints on the operators in Eq. (^) due to measurements of semileptonic 
b — > c decays. They will allow us to bound the coefficient of the operator 0™ R , which contains a right handed charm 
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field and is weakly constrained otherwise. We focus on three types of constraints coming from the ratio of exclusive 
D and D* rates, the polarization in the D* mode, and moments in inclusive spectra. 

We begin with the exclusive case where the B — » Dtv and B — > D*tv rates can be calculated in an expansion in 
Aqcd /n^b,c using heavy quark effective theory. The form factors at zero recoil, where w = v ■ v' = 1 (v and v' are the 
four- velocities of the B and mesons, respectively), have been determined from lattice QCD pc[| . In the SM the 
ratio of rates is independent of V c b, and therefore it provides a good test for non-SM contributions. The presence of 
the new operator, 0™ R , affects the two rates differently. The rates are given by 



dw 487r 3 

(1 _ r*f + (1 - 2wr* + r* 2 ) 

v 1 + w y ' 



dT(B^D*£u) Glm% , 3 >-= - M n2 

; - 1 B r* 3 Vw 2 -l(l + w) 2 



dw 487r 3 



\V cb \ 2 (F D *) 2 , (17) 



where r = mo/ms and r* — mo* /tub- The form factors Tt> & n d J~v* can be decomposed in terms of 6 form factors, 
h+ t - t v,Ai,A2,A3 IM - At leading order in the heavy quark limit !F(w) — J-*(w) = h+ t v,Ai,A 3 — 6( w )j where £,(w) is the 
Isgur-Wise function p2j , while h-^ 2 = 0. Therefore, it is useful to define the following ratios of form factors 

d r \ hv n t \ h A3 +r*h A2 , . 

Ri{w) — - — , R 2 (w) = - , (18) 

which are equal to unity in the heavy quark limit and have been measured experimentally. 

Following the analysis of we can absorb the new physics contributions in the form factors. We obtain 

Ah + = fc(l + r)(l - w)£(w) , Ah- = -k(l - r)(l + u>)f (to) , 

A/iAi = -2fc(l - r*)CH , Aft A2 - -2fe£H , (19) 

A/ia 3 = -2h£(w) , AAv = 2fc(l + r*)CH , 



2vm Bn JC]S R Vtb 



where 

For the new physics contribution we include only the leading term, so we set £(1) = 1. Setting A = 1 TeV, we obtain 



Fd{1) * 


, ^B M (i) 


- 1.01 x 10" 


3 x Re 




* ^(i) 


- 2.02 x 10" 


3 x Re 


fli(l) ? 


a E? M (1) 


+ 6.52 x 10" 


3 x Re 


i? 2 (l) * 


* ^M(l) 


- 2.48 x 10" 


3 x Re 




(21) 

Recent lattice QCD calculations @ give ^|, M (1) = 1.074 ± 0.024 and J^{1) = 0.91 ± 0.04. For Rf M and i?f M we 
use the results of |^4| , scanning over the hadronic parameters that enter. The experimental results are | Vcbl-^ZJ (1) = 
(42.4 ±4.4) x 10~ 3 , \V cb \F D *{l) = (36.2 ±0.6) x 10~ 3 0, i?i(l) = (1.417± 0.075), and R 2 (l) = (0.836± 0.043) ||. 
We set \Vtb\ = 1 and do a combined fit for Cf jR and |T4i>|. We find 

_0.2<^1«M<1.6. (22) 
I Kb | 

We next turn to inclusive B — > XJLv decays, which is also sensitive to the presence of the additional operators. We 
concentrate on the partial branching ratio and moments constructed from the charged lepton energy spectrum (see, 
e.g.,@), 

f dr L E t MrdE t [„ [E e — Mi(E Q )] k dEe 
M (E ) = r B ^dEz, Ml (E Q ) = Je ° ™< , M k (E ) = Je ° 1 Ir 'I ' ■ ( 2 3) 
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FIG. 3: Constraints on 0™ R in the Re(C]f i j) - \V c b\ plane from semileptonic B — > X c £v (solid curves) and B — > D^tv decays 
(dashed curves) and their combination (shaded areas). For each constraint the 68%, 95% and 99% CL regions are shown. 



These are well measured and can be reliably calculated. We use the SM prediction inclu ding l/m 2 and a s corrections 
and compare it in a combined fit with the 20 Babar |27| and a subset |2^] of the 45 Belle [ p9[ measurements, including 
their correlations. The modification of dT/dEi due to the C™ R coupling is given by 

dr NP (5 -> X c iv) = _ G 5 J 2 m e b v 2 Re(C^ R V cb ) ^py 2 (3 - y)(l - y - p) 3 
dy 6\/2ir 3 A 2 (1 - yf 

V2GW b v 4 \CI R \ 2 y 2 (3-y)(l-y-p) 4 
+ 3tt3A 4 (1-y) 3 ' 1 ' 

where y = 2Eg/irib and p — m 2 /m 2 . It is known that the data cannot be fitted well with the OPE truncated at 
l/m 2 . Including the l/m 3 corrections in a more complicated fit would make the agreement with the SM better, and 
therefore our bounds stronger. 

The combined constraints on C™ R and \ V&\ from exclusive and inclusive decays is shown in Fig. [|. The solid curves 
show the constraints from inclusive decays, the dashed curves show the bounds from exclusive semileptonic decays to 
D and D*, and the shaded regions show the combined constraints (the confidence levels are as in Fig. ||). 



D. Exclusive and inclusive b — ► ulv decays 



We now turn to some 3rd — > 1st generation transitions. While the experimental constraints are less precise for these 
than for 3rd — > 2nd generation transitions, the SM also predicts smaller rates, and therefore NP could more effectively 
compete with the SM processes. These constraints are particularly important as they bound the 0[ contributions 
relevant for t — > u decays, which might not be distinguishable at the LHC from t — > c. 

As is the case for 3rd — > 2nd generation transitions, exclusive and inclusive semileptonic b — > u decays can constrain 
the operator 0'^ R in t — > u transition. Similarly to b — > c£D, this operator distorts the lepton energy spectrum, so 
information on the lepton energy moments could constrain it. However, such measurements are not yet available for 
B — > X u lv. Therefore, to distinguish between the SM V u b contribution and C™ R , we use B — * tt£D in addition to the 
inclusive data. 

For exclusive B — > iriP decay, we use for the SM prediction the parameterization of Ref. which relies on 
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FIG. 4: Constraints on 0'™ R in the Re(C2jj) - |Vtii>| plane from B — > X u lv (solid curves) and B — > (dashed curves) and 
their combination (shaded areas). For each constraint the 68%, 95% and 99% CL regions are shown. 

analyticity constraints and lattice QCD calculations of the form factors at large q 2 [}2| . The NP contribution is 



dq 2 24tt 3 1 A 4 

4m B vRe(V tb CZ R V: b ) 
A 2 



[(l-q 2 )fl + (l + f)f-f + ]} 



(25) 



where the f± form factors are functions of the dilepton invariant mass, q 2 , q 2 — q 2 /m 2 B , and we neglected terms 
suppressed by m\jm 2 B . 

For inclusive B — > X u £9 decay wc focus on the measurement utilizing combined cuts (3j| on q 2 and the hadronic 
invariant mass, mx, and compare it with the Belle and Babar measurements [ p4| . Using this determination of V u b 
is particularly simple for our purposes, because in the large q 2 region the mild cut on mx used in the analysis only 
modifies the rate at a subleading level. Working to leading order in the NP contribution, we can neglect the effect of 
the mx cut on the NP and include the NP contribution to the rate via 

dr"(^ x.a» _ ggf |cap ^ _ lf{f + 2) (26) 

Since the interference between the SM and NP is suppressed by m u /mt (see the yfp factor in Eq. (HJ) in the first 
term), there is no dependence on the weak phase of C'^ R in the inclusive decay. Using other determinations of V ub 
would be harder to implement and would not change our results significantly. 

The combined constraint on C' T % a nd \V u b\ from inclusive and exclusive decays is shown in Fig. |^. (This uses the 
lattice QCD input from Fermilab jplfl , and the one using the HPQCD calculation J|2| would also be similar.) 



E. B — * p7 and B —y n + 



A* 



The inclusive B — > Xdj decay has not been measured yet, and there is only limited data on B — > fry. Averaging 
the measurements |35) using the isospin-inspired 4 relation B(B — > pry) = B(B ± — > j o ± 7) = 2(r B ± /rgo) B(B° — > p 0r y), 



4 Isospin is not a symmetry of the electromagnetic interaction. This average relies on the heavy quark limit to argue that the dominant 
isospin violation is Aqcd/"*6 suppressed. With more precise data, using only B° decays will be theoretically cleaner, because annihilation 
is suppressed in the B° compared to the modes. At present, this would double the experimental error, so we include the B^ data. 
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and the PDG value t b ± /t b o = 1.07, we obtain 



B{B -> fry) = (1.26 ±0.23) x 1(T 6 . (27) 



To reduce the sensitivity to form factor models, we normalize this rate to B(B — » K*j) 
[B{B ± -> lf* ± 7) + (t b ± /r B o) 6(5° -» #*° 7 )] /2 = (41.4 ± 1.7) x 10~ 6 , 



BjB -> P7 ) 
6(5 -> if * 7 ) 



/,</ 



2/9 9 



We use £ 7 = 1.2 ± 0.15, where this error estimate accounts for the fact that we consider the rates to be determined 
by Oi-y(mb) alone. The contributions of other operators have larger hadronic uncertainties and are expected to 
partially cancel |}6). If first principles lattice QCD calculations of the form factor become available then one can 
avoid taking the ratio in Eq. (]2|), and directly compare the calculation of B(B — + pj) with data. We obtain the 
following constraints 

- 0.26 < C'£ L < -0.21 or -0.026 < C'£ L < 0.03 , 
-1.2 < C'[\ < -0.9 or -0.11 < C'[ l L < 0.13 , 
-0.7 < C'n L < -0.5 or -0.07 < C'g L < 0.08 , 

-0.1 < C% L < 0.09 or 0.7 < C RL < 0.9 . (29) 
Note that there are no constraints on 0'^ R or 0'^ R , because of their m u /mw suppression. As for B — > A s7 , the two 



solutions in Eq. (29) correspond to the sign ambiguity in interpreting the constraint on |C7 7 | 2 when we assume that 
the NP contributions are real relative to the SM. Had we not made this assumption, the allowed regions would be 
annuli in the complex C- planes. 

The NP operators we consider also contribute to the rare decays Bd, s — > M + M - This is most interesting for 
B c [ — > since one expects that the NP contribution is enhanced compared to the SM by [(w 2 /A 2 )(l/|V t( j|)] 2 , 

which is around 20 for A = 1 TeV. Moreover, 0'^ L contributes at tree level, so its contribution is enhanced by an 
additional factor of (47r/a) 2 . Although this decay mode has not yet been observed and the present upper bound 
B(B — > /i + /i~) < 3 x 10~ 8 |37j is two orders of magnitude above the SM expectation, it still gives a useful constraint 
on 0'£ L . In particular, for A = 1 TeV, we obtain 

- 0.023 < C'£ L < 0.026 . (30) 

The combined constraints from B — * p 7 and B — * on 0'^ L and 0'^ L are shown in Figure |[ The region between 

the dashed lines is beyond the LHC reach, but the LHC will be able to exclude (though perhaps not completely) the 
non-SM region in Fig. g|. In the case of 0'^ L and 0' RL the present data are not strong enough to exclude the non-SM 
region allowed by B — > pj. 



F. AF = 2 transitions 



In this section we present the results of the analysis of the NP effects in AF = 2 processes. Since their contribution 
appear at the same time in B^Bd, B S B S and K°K°, we pe rformed a full fit using the CKMFitter code 0, after 



having suitably modified it to include the results of Sec. Ill C . The code simultaneously fits experimental data for the 



Wolfenstein parameters 7j and rj and for NP (extending earlier studies in AF = 2 processes |38|, |39| ) . The observables 
used here include the Bd and B s mass differences, the time dependent CP asymmetries in B — > J/tJjK, the CP 
asymmetries in B — > irn, pp, p-K, the ratio of \V u b\ and \V c b\ measured in semileptonic B decays, the CP asymmetries 
in B — > DK, the width difference in the B S B S system, Ar s , the semileptonic CP asymmetry in B decays, Asl, and 
the indirect CP violation in K decays, €k- We allowed the NP Wilson coefficients to be complex and performed a 
scan over their phases. Thus, the results in this section are quoted in terms of the absolute values of the Cj and C[. 
Keeping only one operator at a time, we get 

|C£J<0.07, \C h Lh \< 0.014, |C3k|<0.14, 

| C2* L | < 0.11 , \C'£ L \ < 0.018, \C% L \ < 0.26. (31) 

As before, we also performed a combined analysis for the LL operators. This is particularly interesting for 0'£ L and 
0'l L , since until B — > Xdl + t~ data becomes available, only AF = 2 processes are sensitive to the complex phases. In 
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FIG. 5: Constraints from B — > p7 and B — > in the C^i - C^j, plane. The red, green, and blue regions denote 68%, 95%, 

and 99% CL, respectively. The region between the dashed lines is beyond the LHC sensitivity. 



general, allowing for a variation of the phases of C'^ L and C'^ L , a cancellation can occur between the two contributions 
and the above bounds are relaxed. If their absolute values satisfy \C'^ L \ ^0.1 \C'^ L \ then arbitrarily large values of 
the Wilson coefficients is allowed for some values of the phases. This possibility is ruled out when the B — > pj and 
B — > fi + ji~ constraints are included. Indeed, combining A_F = 2 with these measurements, we obtain 



|C^| < 0.26, 



C'l L \ < 0.026. 



(32) 



V. COMBINED CONSTRAINTS AND CONCLUSION 



In this paper, we studied constraints on flavor-changing neutral current top quark decays, t — ► cZ, «Z, cy, 1*7. 
We used an effective field theory in which beyond the SM physics is integrated out. In the theory with unbroken 
electroweak symmetry the leading contributions to such FCNC top decays come from seven dimension-6 operators of 
Eq. (Q). We assumed that the new physics scale, A, is sufficiently above the electroweak scale, v, to expand in v/A 
and neglect higher dimension operators. We find different and sometimes stronger constraints than starting with an 
effective theory which ignores SU (2)^ invariance. 

The 95% CL constraints on the Wilson coefficients of the operators involving 3rd and 2nd generation fields are 
summarized in Table |. We consider one operator at a time, i.e., that there are no cancellations. The top two rows 
show the present direct constraints and the expected LHC bounds. The next three rows show the bounds from B 
physics. In the B — > X s j, X s £ + i~ row the combined bounds from these processes are shown. The two allowed 
regions are obtained neglecting the complex phases of the operators (see Fig. ^| and the discussion in Sec. 
assumption can be relaxed in the future with more detailed data on B ^ X s e + l~. In the AF = 2 row t 
refer to upper bounds on the magnitudes of the Wilson coefficients and are derived allowing the phase to vary. The 
best bound for each operator is listed and then translated to a lower bound on the scale A (in TeV, assuming that the 
C's are unity), and to the maximal t — > cZ and t — > cy branching ratios still allowed by each operator. The last row 
indicates whether a positive LHC signal could be explained by each of the operators alone. In this row, the star in 
"Closed*" for C^ L and C\ L refers to the fact that small values of these Wilson coefficients cannot give an observable 
top FCNC signal, however, there is an allowed region with cancellations between the SM and the NP, which may still 
give a signal. In the same row "Ajar" means that C% L and C b RL cannot yield an LHC signal in t — > cZ but may 
manifest themselves in t — > cy. It is remarkable that the coefficients of several operators are better constrained by B 
physics than by FCNC top decays at the LHC. 



12 





r~1U 


c h LL 


i~1W 

^RL 


^RL 


^LR 


/~tb 
^LR 


^RR 


direct bound 


9.0 


9.0 


6.3 


6.3 


6.3 


6.3 


9.0 


LHC sensitivity 


0.20 


0.20 


0.15 


0.15 


0.15 


0.15 


0.20 


B — * XJ + £~ 


1— 07 0361 


[-0.017, -0.01J 
[-0.005, 0.003] 


1— 09 181 


1— 12 241 


[—14 71 


1— 10 191 




AF = 2 


0.07 


0.014 


0.14 










semileptonic 










[0.3, 1.7] 






best bound 


0.0 ( 


0.014 


O.lo 


U.24 


1 *7 


0.3 


9.0 


A for d = 1 (min) 


3.9 TeV 


8.3 TeV 


2.6 TeV 


2.0 TeV 


0.8 TeV 


0.4 TeV 


0.3 TeV 
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TABLE I: 95% CL constraints on the Wilson coefficients of the operators involving 3rd and 2nd generation fields for A = 1 TeV. 
The top two rows show the present direct constraints and the expected LHC bounds. The second part shows the bounds from 
B physics, which is then translated to a lower bound on the NP scale, A, and to the maximal t —* cZ and t — > cy branching 
ratios each operator could still give rise to (the ATLAS sensitivity with 100 fb" 1 is 5.5 x 10~ J and 1.2 x 10~ J , respectively). 
The last line concludes whether a positive LHC signal could be explained by each of the operators. 
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TABLE II: Constraints on the Wilson coefficients of the operators involving 3rd and 1st generation fields. The entries in the 
table are analogous to Table |. 



Tabic |0] shows the constraints on the operators involving the 3rd and 1st generation quarks. We studied this 
because the LHC may not be able to distinguish between t — > c and t — * u FCNC decays, and these processes are also 
interesting in their own rights. In this case there are two allowed regions of C'™ R from semileptonic decays, as can be 
seen in Fig. ||. The entrie s in the "combined bound" row show the result of the fit to all the B decay data above it, 
as discussed in Sec. LVF. We see from the last row that the LHC window remains open for all of the RR, LR, and 
RL operators, except 0'^ L . 

We conclude from Tables Q and @ that if the LHC sees FCNC t decays then they must have come from LR or RR 
operators, unless there are cancellations. Moreover, if t — > cZ is seen but t — ► cy is not, then only 0^ R could account 
for the data. 

Our analysis used the currently available data and compared it to an estimate of the LHC reach with 100 fb" 1 . 
However, by that time many of the constraints discussed above will improve, and new measurements will become 
available. The direct bounds will be improved by measurements from Run II of the Tevatron in the near future. All 
the B decay data considered in this paper will improve, and the calculations for many of them may become more 
precise. Important ones (in no particular order) are: (i) improved measurements of B — > X s £ + £~ to better constrain 
the magnitudes and especially the phases of C£ L , C^ L , C RL , and C RL ; (ii) measurements of the lepton energy and 
the hadronic mass moments in B — > X u £& to constrain C^; (hi) improvements in B —> pj and measurement of 
B — > Xd'y to reduce the uncertainties of C'j^ L , C'^ L , C RLl and C RL ; (iv) measurement of and B — > Xd£ + t~ to reduce 
the errors and constrain the weak phases of these last four coefficients. Additional information will also come from 



13 



LHCb. For example, the measurement of the CP violating parameter SB s -npcj>, the direct measurements of the CKM 
angle 7, and some of the above rare decays will help improve the constraints. With several of these measurements 
available, one can try to relax the no-cancellation assumption employed throughout our analysis. Note that not 
all NP-sensitive B factory measurements can be connected to FCNC top decays; e.g., the CP asymmetry Sk*j is 
sensitive to right-handed currents in the down sector and cannot receive a sizable enhancement from the operators in 
Eq. (|J). Thus, there are many ways in which there can be interesting interplays between measurements of or bounds 
on FCNC t and b quark decays. 

If an FCNC top decay signal is observed at the LHC, the next question will be how to learn more about the 
underlying physics responsible for it. With a few tens of events one can start to do an angular analysis or study an 
integrated polarization asymmetry |40|j . These could discriminate left-handed or right-handed operators (say O rr or 
0^ L ). Such interactions could arise in models in which the top sector has a large coupling to a new physics sector, 
predominantly through right-handed couplings pd"| . However, a full angular analysis that could also distinguish O rr 
from 0™ R requires large statistics, which is probably beyond the reach of the LHC. 

The observation of FCNC top decays at the LHC would be a clear discovery of new physics, and therefore it would 
be extremely exciting. Our analysis shows that an LHC signal requires A to be less than a few TeV. This generically 
implies the presence of new particles with significant coupling to the top sector. If the new particles are colored, 
we expect that they will be discovered at the LHC. It would be gratifying to decipher the underlying structure of 
new physics from simultaneous information from top and bottom quark decays and direct observations of new heavy 
particles at the LHC. 
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APPENDIX A: ANALYTIC EXPRESSIONS 

We give the form of the operators of Eq. (j^) after electroweak symmetry breaking, keeping only trilinear vertices 
which do not involve the Higgs: 

\f2m 2 w , - . s 2m z m w 
°ll = — WW s l + b L W c L ) H t L Z C L + ... , 

" 92 92 

h 2m z m w , 

O ll = {t L Z c L + b L Zs L ) + ... , 

92 

Orl = mwSLO-^tRW^ + V2m w CL(T^t R (c^Z^u + s w A^ v ) + . . . , 
Orl = V2m w CL^ v t R s w A^ v - — 

°lr = m w b L a^ v c R W~ v + V^mwtLU^CR (c w Z^ v + s w A^) 
- ( s2 

°LR = VZmwtLCT^CR SyjA^ Zp V 

RR = 2 -^^t R Zc R + .... (Al) 

92 

Here s w = sin9 w , c w — cos8 w , and the dots denote hermitian conjugate and the neglected vertices involving Higgs 
and higher number of fields. Throughout this paper the covariant derivative is defined as = 9 M + igA^T a + ig'B^. 
The analytic expressions for the contributions of the operators in Eq. (||) to the top FCNC partial widths are 

T(t - cZ) = (1 - yf{ [\C h LL + CW + \CU 2 ] (1 + 2y) 



16tt A 4 

2g 2 2 cos 2 6 W (2 + y) [\C h LR t<m 2 9 W - C% R \ 2 + \C b RL tan 2 9 W - C£ L ' 2 " 
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6V2 32 sine w tan Q Ws Jy Re [(C^)*(C^ + C£ L ) - (C£ R )*C£ 
672 ff2 cos [{C^ L T{C h LL + C U LL ) - (C^ R )*C U RR ] } 



T(t -> c 7 ) 



2 9 

irmf 



C v r 



CI 



(A2) 
(A3) 



where y — m%/m^. The analogous expressions for t — > u decays are obtained by replacing Ci by C[ above. This 
expression makes it straightforward to relate the Wilson coefficients used in this paper with different notation present 
in the literature, which defines the couplings in the effective Lagrangian after electroweak symmetry breaking. 

Next we present the analytic expression for the Wilson coefficients originating from the operators in Eq. (j^). We 
use the MS scheme and match at the scale /i = mxy- It is easiest to express the results as modifying the Inami- 
Lim functions B , Co, and D /D' Q , coming from box diagrams, Z-penguins, and 7-penguins, respectively. Using the 
standard normalization of the effective Hamiltonian 



V tb V t * s J2dOi, 



(A4) 



CV 7 

C10A — 



'Mx) 



the Wilson coefficients at the matching scale can be written as 

1 

~2 

a [ 1 

2n \_ sin 2 9w 
a 1 



B (x) 



1 



sin 6w 



4 )Co(x)-D (x) 



2tt 



sin &w 



[B (x)-C {x)], 



(A5) 



where x = m^fw^. In the SM, we have the well-known expressions [ fL2[ 



B (x) 
C (x) 
D {x) 
D' (x) 
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In a; , 
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(A6) 



The contributions of the 0^ L , 0^ R , and O h LR operators introduced in Eq. 
terms to Eq. ( A6 ) 
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Note that the contribution of 0\ L to A Co occurs at tree level, as indicated by its 1/0:2 enhancement in Eq. flA8|), so 
0\ L gives tree- level contributions to Cgy and C\qa- Nevertheless, we shall not include the matrix element of Oj jL to 
one higher order in 0:2, in analogy with the conventional approach in which the NNLL calculation of B — > X s £ + £~ 
does not include the 0(a 2 ) matrix element of Ogv- 

Finally we calcuate the AF = 2 contributions due to C£ L and C™ R . The shift in the SM contributions read 
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The Kij are zero except for K K mixing, where they are given by 
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